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ABSTRACT 


Solar cycle variations of the proton population of the inner radia- 
tion belt are calculated in this report. The analysis includes proton 
energies from 10 to 700 Mev, field lines from 1.2 to 1.6 earth radii, and 
appropriate values of field strength. An averaged atmospheric model 
is constructed in terms of R (flux density in gauss) and L (earth radii) 
coordinates which represents the average proton number density that 
trapped particles would encounter. The averaging process includes 
diurnal, longitudinal, north-south, and "bounce" averaging. The time 
dependence of the atmospheric model is constructed for the current 
cycle using averaged information of the 10. 7-cm solar flux of the recent 
past. Proton flux spectra for solar minimum and solar maximum are 
calculated from the particle conservation equation. Transient proton 
spectra and time histories are presented showing the dynamical be- 
havior of trapped protons when influenced by the fluctuating atmosphere 
and a variable source. The results indicate that in the region of B-L 
space corresponding to minimum altitudes of 300 to 800 km, the proton 
flux is changed appreciably from solar maximum to solar minimum 
where solar minimum is the larger. Calculations indicate a change in 
the nature of the proton spectra where the spectra become peaked at 
some energy level depending on the time in the cycle and position in 
space. Comparisons are made with available data. 
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SOLAR CYCLE EFFECTS ON INNER ZONE PROTONS 


by 

R. C. Blanchard 
and W. N. Hess 
Goddard Space Flight Center 


INTRODUCTION 

The purpose of this paper is to calculate the expected changes in inner-zone proton densities 
with time in the solar cycle. Freden and White identified the penetrating component of the inner 
belt as due to energetic protons (Reference 1). Subsequent experimental work by Freden and White 
(References 2 and 3), and Armstrong, Harrison, Heckman, and Rosen (Reference 4) has extended 
this finding and we now have a well-established proton energy spectrum at earth radii, L, approx- 
imately equal to 1.4, and flux density, B, approximately equal to 0.20 for a time near solar maxi- 
mum (Figure 1). Their measurements have established the energy spectrum of the protons for 
E > 20 Mev. The analysis of this spectrum (Reference 3) by an albedo neutron decay source 
mechanism suggests that this source adequately describes the production of protons of the inner 
zone. 

One of the features of this component of the belt is its time consistency. Measurements after 
a solar flare by Freden and White (Reference 3) indicated little variation in proton flux. Yoshida, 
Ludwig, and Van Allen (Reference 5) have shown from Explorer I data that the change in the proton 
population is less than a factor of 2 for a two- month period. Data from Explorer IV (Reference 6) 
indicated the same results for a three- month period in 1958. Data from Explorer VII (Reference 
7) collected for a 15-month period indicated a net change by a factor of about 4 or 5 at the outer 
part of the inner zone to about 2 or 3 at the inner part of the inner zone. Although the net change 
in Explorer VII data is not pronounced, it does show a steady increase in the counting rate with 
some superimposed variations. Hess (Reference 8) has shown that slow changes in proton popula- 
tion were expected due to solar- cycle effects on the galactic cosmic ray flux and, more importantly, 
on the fluctuating atmosphere due to solar exospheric heating. It turns out that these two factors 
are resonant, producing a larger proton flux at solar minimum than at solar maximum. Hess (Ref- 
erence 9) estimated this factor to be 10 or more over a solar cycle. It is believed that some of the 
steady increase in the counting rate from Explorer VH at a time approaching solar minimum could 
be due to the solar-cycle effect. This report will examine the solar- cycle effect on inner- zone 
protons to determine its characteristics in detail. 
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CONTINUITY EQUATION 

The study of the variation in the proton 
population requires knowledge of the conditions 
which govern the source and loss of protons in 
a given region. The continuity equation for 
inner- zone protons is usually written as (Ref- 
erences 2, 10 and 11) 

dN(E) d / dE\ 

“dT = S( £ ) ' L ( E ) + dE (Nv * (1) 

where S(E) is the albedo neutron decay source 
of protons being supplied to the region; L(E) is 
the nuclear interaction loss of protons in the 
region; and the last term includes the effect of 
proton energy loss due to near coulomb colli- 
sions with bound atmospheric atoms. This 
equation approximately represents the conser- 
vation condition of trapped protons, which yields 
the rate of change of the proton number density, 
while integration of this equation gives the total 
proton number density at a given energy level, 
E, in protons/cm 3 -Mev. 

The loss term L used in this study is basi- 
cally that used by Freden and White (Reference 
2) with some modifications. The term takes into 
account catastrophic nuclear collisions of trap- 
ped protons with the atmospheric atoms. This 
term is written as 



Figure 1— The experimentally measured inner-zone pro- 
ton energy spectrum (Freden and White, Reference 3) for 
L ~ 1.4, B ~ 0.20 at time near solar maximum. 


5 

L = Nv^]nicr(j) , (2) 

j = l 

where j refers to the atmospheric components considered, which are 0 2 , O, N 2 , He, and H, n j is 
the average atmospheric density of the j th constituent and cr( j ) is the inelastic geometric cross 
section of the j th constituent. By defining 

5 

j = 2 
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Equation 2 can be written more conveniently as 


L = NvS 


(4) 


The values of or for oxygen and helium used in this study are 

cr(O) - 0.36 barns 


and 


cr(He) 


0. 143 barns 


For simplicity of calculation, it is assumed that 
the nitrogen interaction cross section is equal 
to that of oxygen and the hydrogen contribution 
is negligible. Thus, the calculation of 2 is made 
by the relation: 


2 $ 


l (He) 


2 cr(He) + 


n<°> + 2n ( ° 2> 


2 s (N2 > 


CT (0)- (5) 


The pulsation of the atmosphere causes s 
to be a function of time in the solar cycle as 
well as position. Figure 2 shows a time history 
of log e 2 for various values of B at L = 1.25. It 
can be seen that as B increases, the variation 
from solar maximum to solar minimum de- 
creases. The fluctuating atmosphere is not as 
pronounced at large flux densities which cor- 
respond to low altitudes. 



Figure 2—A time history of the "effective" cross section 
of the atmosphere, 2 as a function of B at L = 1.25 e.r. 


The source term S(E) used in this study is 
from albedo neutron decays. The neutrons are 
produced from cosmic ray protons colliding 
with atmospheric oxygen and nitrogen. The 

produced neutrons diffuse out from the atmosphere and subsequently decay into protons after ap- 
proximately 10 3 seconds. The neutrons which escape from the atmosphere and then decay 
will be injected into the belt. The form of the proton source to be used is essentially the 


3 



same as that used by Hess (Reference 11) with some minor modifications. Specifically, these 
are the addition of the solar cycle time dependence and the transformation into B-L space. 
This leaves the source as: 


0 . 8E~ 2 $ 
fiCyrh 2 cos 4 k 


(6) 


where 

v = neutron velocity (cm/sec), 

P = v/c , 

C = speed of flight, 2.9979 X 10 10 cm/sec, 

> = (l-/3 2 )~ 1/2 , 

r = neutron mean life, 10 3 sec, and 
0 . 8E~ 2 = calculated neutron flux (neutrons/cm 2 -Mev~sec), 


and where the magnetic latitude dependence, is replaced by B through the dipole magnetic field 
equation. This cos -4 k term is left in the equation due to the inability of the field equation to be 
solved readily in closed form. The nondimensional parameter o is the relative inner belt source 
strength. That is, 


4 >( t, E) 

*(V E ) ’ (7) 


-e* 



Figure 3— A time history of the relative inner-belt 
source strength. 


where the nomalizing factor $(t 0 E) = 0.8E" 2 
neutrons/cm 2 -Mev-sec. The relative inner- 
belt source strength is a function of time in the 
solar cycle. Figure 3 is a plot of a current 
estimate* (Reference 12) of the change of $ for 
the current solar cycle. Three points of inter- 
est in this graph are: first, the net change over 
the cycle is 25%; second, at solar maximum the 
production of neutrons is smallest due to ex- 
clusion of some neutron- producing galactic 
particles by the sun r s increased activity; and 
third, $ is not symmetrical. 

It is assumed that the source, s, produces 
protons such that all produced protons have 


♦McDonald, F. B., private communication. 
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velocity orientations perpendicular to the field at the magnetic mirror latitude, This assumption 
is made instead of adding the contributions of protons at other positions along a field line which 
have the necessary mirror point conditions (that is, the proper pitch angle). The assumption of an 
injection coefficient of 1 will probably not affect the overall result of the study, which is the rela- 
tive change from solar minimum to solar maximum. 

Equations 6 and 4 represent the form of the components of Equation 1. Approximate values 
of ft and y were generated by Freden and White (Reference 2). They are: 


ft = 0. 0484E* 477 

7 = 0. 930E' 032 


for 10 <E< 80 Mev 


ft = 0.0896E* 344 

7 = 0. 428E* 205 


for 80 < E < 700 Mev 


Utilizing these approximations and Equations 4 and 6, Equation 1 can be arranged in the following 
form: 


dN 

dt 


a n ® 


L 2 E 0 cos 4 A. 


E b ’ 


b 

, NE 2 


A (A\ 

dE \dX/ 


, E b2 N 2 


where, if the energy is 10 < E < 80 Mev 




a o = 

2.694 x 

10“ 13 

b o 

= 2.509 

a l = 

3.463 x 

10 8 


= 0.523 

a 2 = 

7. 255 x 

10 s 

b 2 

= 0.477 

and if the energy is 80 <E <700 Mev 





a o = 

3.479 x 

10“ 13 

b o 

= 2. 549 

a l = 

4.617 x 

10 8 

b i 

= 0.656 

3 2 = 

1.343 x 

10 9 

b 2 

= 0.344 


The calculation of the energy loss due to slowing down is approximately evaluated by 



( 8 ) 


(9) 
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where the scale factor, R, is the ratio of the equivalent oxygen number density of the atmosphere 
to the oxygen number density at NTP conditions. That is, 

1.75n^ N ^ + 5<°> + 0.25n< He > + 2n^° 2 ^ + 0. 125S< H) 

R 2.69 x 10 19 


Figure 4 is a time history of R at L = 1.25 for various B’s. Again the solar-cycle variation be- 
comes apparent. The values of (dE/dX) NTP comes from experimental data published by Aron, Hoff- 
man, and Williams (Reference 13) for protons penetrating an oxygen absorber. Figure 5 is the 
energy spectrum of this data for the energies of interest. In a similar fashion the energy slope of 
the energy loss is 


_d_ 

dE 




(10) 



Figure 4— A time history of the atmosphere scale factor, 
R, as a function of B at L = 1.25 e.r. 


where the slope of dE/dX has been calculated 
from Figure 5 and represented in Figure 6. 
Using Equations 9 and 11, Equation 8 can be 



Figure 5— The proton energy loss spectrum for 
an oxygen target. 
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written as 


dN 

dt 


a Q 0 


L 2 E 0 cos 4 K 


1 m /dE_\ 
b i \dX / N TP 


- a 2 NE b2 R 


_d^ 

dE 



N2 


, (ID 


where the constants are those defined previ- 
ously. A Fortran program has been written to 
numerically integrate Equation 11 where the 
functional dependence of the proton density rate 
is f (N, E, L, B, t ) . 

A difficulty arises when evaluating the 
transient proton number density. As is evident 
from Equation 11, an initial value of N is neces- 
sary in order to integrate the equation. At time 
t 0 there are no trapped protons and the problem 
is started at solar minimum. This is done in 
order to investigate the most rapid build-up, 
since at solar minimum the injection of protons 
is the largest while the removal process is the 
smallest. 



Figure 6 — The slope of proton energy loss versus energy 
for an oxygen target. 


Following the work done earlier (References 2, 10 and 11), the steady- state proton population 
is found by assuming the rate of change of the proton number density to be zero. The result- 
ing solution of Equation 11 becomes: 



+ a 2 E 


(§) 


+ a 2 E 2 £ > cos 4 K n 


(12) 


where the coefficients are those defined previously. As can be seen from the functional dependence 
of N, a time must be chosen in order to evaluate the steady- state spectrum for a given position. 

Two times are chosen: maximum solar and solar minimum. This is done in order to investigate 
the maximum change in proton number density at the two source and loss extremes. 


AVERAGE ATMOSPHERE IN B-L SPACE 

A major part of this problem was to determine the average atmospheric density used in calcu- 
lating the scale factor, R, which is contained in the loss- rate terms of trapped protons, and 2, which 
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is contained in the nuclear interaction loss term. The calculation of these quantities has been 
done such that over a long sampling period, on the order of months, the values obtained would be 
representative of what a trapped particle would "see” while moving about the earth's magnetic 
field. Since the particle's three basic motions (that is, gyration around the field line, cyclic motion 
parallel to the field line and drift motion in longitude) are separated by large time differences (Ref- 
erence 14), it is possible to consider them separately in the averaging process. 

The model of the upper atmosphere used in the study is the one developed by Harris and 
Priester (Reference 15). This model gives the time dependences of the atmosphere for both solar 
cycle and diurnal variations. The atomic densities, n., of five atmospheric constituents (i.e., 0 2 , 
O, He, N 2 , H) are given in the form 


n. - n. (h, t, S) , 

where h is the altitude above the earth, t is local time given in hourly increments and S is a model 
parameter related to, but not the same as the intensity of the 10.7 cm solar flux F in watts/ 
(m 2 - cps x 10 -22 ). Harris and Priester's atmospheric model agrees well with several measure- 
ments of density by satellite drag (References 16 and 17) and has been checked by preliminary 
density data from the Explorer XVH satellite (References 18, 19 and 20). To date, this model is 
the most complete description of time dependence of atmospheric densities available and agrees 
well with current experimental data. 

The first step in arriving at the average models is to calculate the diurnal average number 
density of each of the five S models presented by Harris and Priester. The sum, 

2 4 

Hi < s > = h J2 n * <s) . (i3) 

i = 1 

is computed, where the five atmospheric elements are j = He, O, 0 2 , N 2 , H and where i refers 
to the hourly value of the density. The diurnal average is taken, since protons considered in this 
study drifting in longitude around the earth have periods of revolution on the order of 1 to 30 
minutes (Reference 14). Thus, over a period of 24 hours, the daily proton population variation will 
tend to be averaged out, or at least to be a second- order effect compared to the solar- cycle ex- 
pansion of the atmosphere which has a period of approximately 11 years. In essence, this step 
eliminates short-term time effects. 

The second step in constructing the solar-cycle average atmosphere is to consider the longi- 
tudinal drift of protons along a B-L contour. As the particles drift around the earth, their mirror 
points encounter different altitudes, and, consequently, different atmospheric densities. To ac- 
count for this, the calculation of the longitudinal average density, n k J , is done by the sum 


n 


j 

k 



i = 1 


n i j (0) 


(14) 
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where K refers to either Northern or Southern Hemisphere and where n. j (0) is the diurnal average 
number density of the j th constituent at longitude <f>. The factor 1/35 is used since equal incre- 
ments of 10° in longitude were used to evaluate the summation. To obtain values of n j as a function 
of longitude in each hemisphere, B-L contours (Reference 6) were generated (Reference 22) as a 
function of altitude and longitude. Figures 7 and 8 show these contours for a typical line of force, 

L = 1.25. Values of atmospheric number density as a function of longitude for each 10° were ob- 
tained by interpolating the results of Equation 13 for altitude which from Figures 7 and 8 corresponds 
to a given longitude for a given B and L. Upon evaluating Equation 14, the resulting functional de- 
pendence of the longitudinal-averaged atmosphere is 

n K J = (B, L, S) 

Both hemispheres are then averaged. As can be seen from Figures 7 and 8, the Southern Hemis- 
phere greatly influences the resultant atmosphere since it dips lowest in altitude due to the nature 
of the earth's magnetic field. Figure 9 is a plot of minimum altitudes for a B-L contour in the 
Southern Hemisphere. Plotting minimum altitude contours in B-L space indicates the region in 



Figure 7— B (gauss) contours for L = 1.25 for the Southern 
Hemisphere as a function of altitude and geographic 
longitude. 



Figure 8— B (gauss) contours at L = 1.25 e.r. for the 
Northern Hemisphere of altitude and geographic 
longitude. 
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L ( earth radii ) 

Figure 9— Minimum altitude contours in B, L space 
for the Southern Hemisphere. 



terms of B and L in which the atmosphere is ap- 
preciable as far as trapped particles are con- 
cerned. The assumption is made that protons 
drift with uniform velocity in longitude. This 
eliminates weighing the average process over 
longitude. Recently, Walt* indicated this is not 
the case. Thus, some error is introduced by 
not appropriately weighing the longitudinal 
averaging; probably about 20%. 

The last averaging step is to adjust the 
longitudinal average density n j (B, L, S) due to 
the protons' north-south mirroring motion. 
This process reduces the magnitude of the 
tables due to motion of a proton from a mirror 
point at low altitude and high density moving 
towards the equator at a high altitude and low 
density. Since the protons live much longer 
than a bounce period, and since no change in 
direction of the proton is considered during 
slowing, it is permissible to average the den- 
sity in this fashion. As one might expect, the 
account by which n J is reduced depends upon 
the arc distance away from the equator; as the 
arc distances approaches zero, the correction 
becomes unity. Following Ray (Reference 23), 
the bounce average of the number density is de- 
fined as 

n(B, L, S) ds 



Figure 10-Schematic of a trapped particle's north-south 

motion where ds is the element of arc along the par- where ds is “ element of arc along the parti- 
tide's helical trajectory, dl is the element of arc along cle's orbit. Figure 10 is a schematic of a 
the field line. trapped particle's north-south motion. The earth 

is assumed to be a dipole for the calculation of 
this averaging step. Due to symmetry, the integral need only be evaluated over a quarter of a 
complete oscillation. The procedure adapted for the calculation of Equation 15 is to project the 
element of arc onto the field line. This is done for convenience since the atmosphere is given in 
terms of field lines, L . The element of length along the particle's orbit is related to an element 


*WaIt, M., private communication. 
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of length along a field line by 


vdt _ dl 

c * s ” cos a ~ cos a ’ 

where a is the pitch angle (i.e., the angle between the velocity v and the magnetic field b). By 
combining 


dl 2 - dr 2 + r 2 dA 2 


along with the field equations 


B - — “ (l + 3 sin 2 a) 1/2 

r = L cos 2 A , 


and a relationship of the conservation of the magnetic moment, which is 


where 


constant 


reduces Equation 16 to 


rk 

I. * 
f 


J (A) A (A) d k 


A (A) dl 


A (A) - cos 4 A. 


4 ”3 cos 2 A 

i/f - 3 cos 2 k cos 6 A - cos 6 A V4 “ 3 cos 2 A 

f m m f 


1/2 


(16) 


(17) 


The subscript m corresponds to the particle's mirror point. The weighing factor, A(A), appears in 
the averaging equations due to the fact that particles spiral about the field line in such a fashion 
as to stay longer at some latitudes, namely, near mirror latitudes, A m . Figure 11 is a plot of A(A) 
versus A for different mirror latitudes, where the mirror latitudes (A m ) are the asymptotes of each 
curve. A dipping of the curves occurs at A > 48°for large values of A m . This phenomena occurs 
principally because the lines of force become relatively steep for large mirror altitudes. For an 
equal AA, it turns out that the particle will spend less time at latitudes >48°, and not near the 
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Figure 1 1 — The weighing factor, A(k) f versus latitude 
for various mirror latitudes, k m/ where the mirror lati- 
tudes (kj are the asymptotes of each curve. 


mirror point, than at latitudes less than 48 °. To 
elaborate, if one assumes that the mirror lati- 
tude is very nearly equal to 90°, then from 
Equation 17 

. cos \ V4 _ 3 cos 2 X. 

A (M * 17414 * 

To find the point of inflection of A(\), the deriv- 
ative with respect to k is set to zero. Solving 
the resulting equation, it turns out that 


Checking Figure 11, one sees that in this region 
the curve begins to dip. To find the time spent 
per unit path length consider 


where V (( is the component of the velocity along 
the field line. By comparing Equations 15 and 
16, the time spent per unit path length in terms 
of the weighting factor is 

dk 1 ds _ v cos a 

dt " A(\) dt " A(\) 

For particles mirroring at large latitudes, the 
length of path per unit time at the equator is 

AA. % v 


Whereas at k > 48° where the pitch angle is still very near zero, the time spent per path length is 
greater than v since A (k >48°) is less than one. 

The integrals of Equation 16 are numerically evaluated on a digital computer using Simpson r s 
technique. At mirror latitudes, k m , the expression A (k m ) is undefined. In order to overcome this 
difficulty, Equation 16 is numerically integrated to A. m - e where e is made arbitrarily small such 
that the values of the integrals do not change appreciably. 

The last step is to eliminate the model parameter, S, for time in the solar cycle. A current 
estimate of the present solar cycle is given by Figure 12. This has been constructed from infor- 
mation obtained from Harris and Priester (Reference 15). An extrapolation beyond 1961 has been 
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made to order to estimate the flux at the pres- 
ent time which is near solar minimum. The 
constructed mean solar cycle with epoch at 
January, 1954 is assumed to be 11 years dura- 
tion. The unsymmetrical nature of the cycle 
will produce interesting effects on the proton 
population, as will be seen in the following dis- 
cussion of the results. 

RESULTS 

Equation 11 is numerically integrated on a 
digital computer using a fourth-order Runge- 
Kutta technique with the initial conditions of 
N(E) = 0 and t = solar minimum. The termi- 
nation criterion is that a given cycle proton 
density is the same as the preceding cycle pro- 
ton density within an accuracy tolerance of less 
than 1%. Figure 13 is a plot of time to reach 
the termination criterion in terms of solar 
cycles (that is, 11 years) versus energy for a 
typical line of force, L = 1.25. As can be seen, 
the time required to build the steady- state con- 
dition is rather long for low values of field 
strength. A build-up time of 10 solar cycles 
corresponds roughly to 700 km in minimum 
altitude. Below a B ^ 0.22 gauss, which cor- 
responds roughly to a minimum altitude of 250 
km, the test for steady- state is difficult to de- 
tect due to the very low proton density and the 
almost constant atmospheric density. 

Table 1 gives the figure numbers of the B- 
L points considered for this report. The first 
figure number refers to the plot of the proton 
flux time history while the second refers to the 
proton flux spectrum plot. 


5 
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Figure 12— A time history of the constructed mean solar 
cycle variation of the Harris and Priester model param- 
eter, S, with reference time \ Q of June, 1954. 



igure 13 — The time required in terms of solar cycles to 
uild steady-state conditions versus energy as a function 
f B for L = 1.25 e.r. 


Figures 14 through 40 are time histories of the proton flux for lines of force, ranging from 
L = 1.188 to 1.600 and various minimum altitudes, after the aforementioned steady-state condition 
has been satisfied. It is seen from these graphs that higher energy protons are not extremely af- 
fected by the solar-cycle variations, whereas the lower energy protons vary considerably with the 
cycle. Due to the "sluggishness” of the higher energy protons, the phase between proton flux 
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Table 1 


Summary of B-L Points. 


Figure Number 

L (e.r.) 

B (gauss) 

K (degs.) 

Kin (kms.) 

Time History 

Energy Spectra 

14 

41 

1.1880 

.1864 

.954 

710 

15 

42 


.1884 

2.954 

650 

16 

43 


.1930 

4.954 

580 

17 

44 


.2075 

8.954 

400 

18 

45 


.2189 

10.954 

290 

19 

46 

1.2500 

.1677 

5.945 

900 

20 

47 


.1741 

7.945 

810 

21 

48 


.1830 

9.945 

710 

22 

49 


.1936 

11.945 

580 

23 

50 


.2070 

13.945 

410 

24 

51 


.2240 

17.945 

210 

25 

52 

1.3419 

.1674 

13.933 

900 

26 

53 


.1809 

15.933 

790 

27 

54 


.1974 

17.933 

550 

28 

55 


.2173 

19.933 

350 

29 

r 

56 

1.4000 

.1661 

16.924 

1050 

30 

57 


.1821 

18.924 

860 

31 

58 


.2014 

20.924 

600 

32 

59 


.2250 

22.924 

330 

33 

60 

1.5000 

.1727 

21.908 

1050 

34 

61 


.1936 

23.908 

750 

35 

62 


.2190 

25.908 

450 

36 

63 


.2501 

27.908 

130 

37 

64 

1.6000 

.1694 

24.905 

1100 

38 

65 


.1926 

26.905 

800 

39 

66 


.2208 

28.905 

480 

40 

67 


.2557 

30.905 

150 
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M ( protons /cm 2 - sec - Mev ) *n M ( protons/cm 2 - sec - Mev ) 



TIME (years) 

gure 26— Proton time histories for L = 1.34, 
B = 0.181, h min = 790. 



Figure 27 — Proton time histories for L = 1.34, 
B = 0.197, h min = 550. 




TIME (years) 


Figure 28— Proton time histories for L = 1.34, 


B = 0.217, h min = 350. 


TIME (years) 

Figure 29 — Proton time histories for L = 1.400, 
B = 0.1661, h min = 1050. 
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Figure 38— Proton time histories for L = 1.6, 
B = .1926, h . = 800. 

• min 


Figure 39 — Proton time histories for L = 1.6, 
B = .2208, h . = 480. 

min 



Figure 40— Proton time histories for L = 1.6, 
B = -2557, h min = 150. 


minimum and solar minimum (i.e., t = 4 years on the 
time scale) is shifted. As can be seen, the amount of 
shift depends on the proton's energy and position in space. 
In some instances this shifting of phase with respect to 
the solar cycle permits lower energy protons to be less 
abundant than higher energy protons, a result not readily 
apparent from the steady- state solution (that is, Equa- 
tion 12). 

Figures 41 through 67 are proton-flux energy spectra 
for different lines of force and various minimum altitudes 
at different times in the solar cycle. The dotted lines are 
the solution of Equation 12 at the two time extremes, 
solar minimum and solar maximum. The aforementioned 
solar- cycle effect of changing the nature of the energy 
spectrum is more apparent from these figures. It is in- 
teresting to note that as the minimum altitude for these 
curves approaches roughly 800 km, the solar cycle varia- 
tion is almost negligible. It is evident now that a more ef- 
ficient way to calculate the proton density in the upper 
region of space where h min > 800 km is to average the 
number density over the solar cycle and apply the 
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Figure 41— Proton energy spectra at 
different times in the solar cycle for 
L = 1.188, B = .1864, h min = 710. 



ENERGY (Mev) 



Figure 42— Proton energy spectra at 
different times in the solar cycle for 
L = 1.188, B = .1884, h min = 650. 




Figure 43— Proton energy spectra at 
different times in the solar cycle for 
L = 1.188, B = .193, h min = 580. 



Figure 44— Proton energy spectra at 
different times in the solar cycle for 
L = 1.188, B = .2075, h min = 400. 


Figure 45— Proton energy spectra at 
different times in the solar cycle for 
L = 1.1880, B = 0.2189, h mjn = 290. 


Figure 46— Proton energy spectra at 
different times in the solar cycle for 
L = 1.25, B = .1677, h . = 900. 

• 'min 


The curves are labeled by the time in years from solar minimum. 
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Figure 47— Proton energy spectra at 
different times in the solar cycle for 
L = 1.25, B = .1741,h mjn =810. 




Figure 48— Proton energy spectra at 
different times in the solar cycle for 
L = 1.25, B = .183,h min = 710. 




Figure 49— Proton energy spectra at 
different times in the solar cycle for 
L = 1.25, B = .1936, h min = 580. 



Figure 50 — Proton energy spectra at 
different times in the solar cycle for 
L = 1.25, B = .207, h min = 410. 

The curves 


Figure 51 — Proton energy spectra at 
different times- in the solar cycle for 
L = 1.25, B = .2240, h mjn = 210. 

re labeled by the time in years from sol 


Figure 52 — Proton energy spectra at 
different times in the solar cycle for 
L = 1.3419, B = .1674, h min = 900. 

r minimum. 
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Figure 53— Proton energy spectra at 
different times in the solar cycle for 

L = 1 .3419, B = .1809, h mjn = 790. 




Figure 54— Proton energy spectra at Figure 55— Proton energy spectra at 
different times in the solar cycle for different times in the solar cycle for 
L = 1.3419, B = .1974, h min = 550. L = 1.3419, B = .2173, h min = 350. 



Figure 56— Proton energy spectra at 
different times in the solar cycle for 
L = 1.400, B = .1661 / h |njn = 1050. 



Figure 57— Proton energy spectra at 
different times in the solar cycle for 
L = 1.400, B = .182, h min = 860. 



Figure 58— Proton energy spectra at 
different times in the solar cycle for 
L = 1.400, B = .2014, h min = 600. 


The curves are labeled by the time in years from solar minimum. 
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Figure 59— Proton energy spectra at 
different times in the solar cycle for 

L = 1.40, B = .225, h min = 330. 



Figure 60 — Proton energy spectra at 
different times in the solar cycle for 
L = 1.5, B = .1727, h min = 1050. 



Figure 61 — Proton energy spectra at 
different times in the solar cycle for 

L = 1.5, B = .1935, h min = 750. 



Figure 62 — Proton energy spectra at 
different times in the solar cycle for 
L = 1.5, B = .2190, h . = 450. 

' mm 



Figure 63— Proton energy spectra at 
different times in the solar cycle for 
L = 1.5, B = 0.2501, h min = 130. 
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Figure 64— Proton energy spectra at 
different times In the solar cycle for 
L = 1-6, B = .1694, h mjn = 1100. 


The curves are labeled by the time in years from solar minimum. 
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Figure 65— Proton energy spectra at Figure 66— Proton energy spectra at Figure 67— Proton energy spectra at 

different times in the solar cycle for different times in the solar cycle for different times in the solar cycle for 

L = 1.6, B = .1926, h. = 800. L = 1.6, B = .2208, h . = 480. L = 1.6, B = .2557, h . = 150. 

The curves are labeled by the time in years from solar minimum. 



Figure 68— The amplitude of the proton population 
changes during the solar cycle for L = 1.25. 


steady- state Equation 12 by averaging over the 
solar cycle. From Figure 13 it is apparent that 
this procedure would eliminate a great deal of 
numerical integration. A general procedure has 
been adopted to utilize this fact by calculating the 
ratio of proton flux at solar minimum to proton 
flux at solar maximum. Figures 68 and 69 are 
contours of this ratio plotted in B versus energy 
for various L lines. When the ratio, which depends 
on the proton’s energy and position in space, is in 
the vicinity of one, the steady- state proton calcu- 
lation with the solar-cycle averaged atmosphere 
will suffice. 


It is seen from Figures 41 through 67 that the 
change in proton density from solar minimum to 
solar maximum is appreciable. The largest variations are found in the lower part of the energy 
spectrum. We may have overestimated the amplitude in this present calculation by as much as a 
factor of 2 because our solar cycle model uses S = 70 as a solar minimum while recent data of 
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I 

| Harris and Priester shows that S = 100 is 

I probably more appropriate. Changing the shape 

I of the mean solar cycle (Figure 12) will also 

| affect the results of the calculation. It is ex- 

| pected that a more symmetrical cycle would 

I affect the shape of the energy spectrums pre- 

| sented in Figures 41 through 67. 

i A comparison of past experimental results 

| and the calculations of this paper is shown in 

Table 2. In discussing these experiments we 
| can consider three types of information: 

? 1. Data obtained from individual experi- 

ments which run long enough (a year or 


Table 2 

Comparison of Experimental Results of Inner Zone Proton 
Measurements with the Calculations in Present Paper. 



Experiment 

Date of 
Experiment 

Time in 
Solar Angle 

L 

B 

H Min. 

Results of Experiment 

Pi edictions from Present Paper 

(a) 

Naugle and Kmffen 

Sept. 1960 

6.7 

(a) 1.54 

.209 

610 km 

Spectrum flat from 20 to 60 MeV 

Spectrum should fall about a factor of 3 in this 
energy range (Fig. 28, curve 7) 


(1963) — Emulsions 
on NERV 






J n (31) at 6.7 

J„ (31) at 4.1 (Exp. IV) k 

J 0 (31) at 6.7 

J 0 (31) at 4 = 13 (Flg - 28) 





(b) 1.47 

.223 

4 00 km 

Spectrum flat from 20 to 100 Mev 

Spectrum should fall about a factor of 4 in this 








J„ (31) at 6 7 

energy range (Fig. 27, curve 7) 








J n (31) at 4 1 (Exp. IV) 2 5 

J 0 (31) at 6.7 

J 0 (31) at 4 ■ 3 5 (Flg ’ 27} 

(b) 

Freden and White 

1960 

6.5 

1 40 

202 

600 km 

Spectrum peaks at 40 Mev 

No spectral peak at t - 6.5 (Fig. 28) 


(1962) - Emulsions 
on Atlas 






J n (20) 

J 0 (40) 75 

J 0 (20) 

J o (40) 2 0 J o (20) 

but at t = 5 do have spectral peak j - 0.9 

(c) 

Pizzella, Mcllwain 
and Van Allen — 

Oct. 1959 
to 

5.8 

to 

1.40 

22 

390 km 

J 0 (E 18) at 7.0 
J 0 (E 18) at 6.0 2 3 

J 0 (E 25) at 7.0 

J 0 (E = 25) at 6.0 = 2 0 (Flg ‘ 24) 


G M counter on 
Exp. VII 

Dec. 1960 

7.0 






(d) 

Filz and Yagoda — 
Emulsions on 

Dec. 1960 
to 

7.0 

to 

varying 

varying 

400 km 

J n (E 55) at 8.5 
J n (E 55) at 7.0 20 

J 0 (E - 75) at 8.5 

J„ (E = 75) at 7.0 23 (Flg ‘ 24) 


Discoverer 

June 1962 

8 5 






(e) 

Heckman and Nakano - 
Emulsions on Low 

Sept. 1962 
to 

8.7 

to 

varying 

varying 

400 km 

J 0 (65) at 9.5 

J 0 (65) at 9T) 10 * 10 

J 0 (65) at 9.5 
J 0 (65) at 9."0 = 12 


altitude polar 

Sept. 1963 

9.7 




J D (65) at 9.5 

J 0 (65) at 9.5 


satellites 






J 0 (65) at 4.1 (Exp. IV) 3 

J 0 (65) at 4.1 = 12 

(0 

Freden and 

Sept, and 

8.8 

1 30 

.200 

530 km 

J 0 (?) 



Paulikos — Solid 

Oct. 1962 





J 0 (85) 13 



state counters on 
Polar satellites 



1.47 

.223 

400 km 

J ° (7) o, 

J 0 (85) 23 

J 0 (10) 

j c 155) 11 (Fig * 27) 





1.54 

.209 

610 km 

J 0 (?) 

Jo (85) ' 43 

J 0 (10) 

Jol^j 24 (Fig - 28) 





low 

high 


J (E • 31) at 8.8 _ 

J (E > 31) at 4.1 (Exo. IV) 3 

J (E = 50) at 8.8 
J (E = 50) at 4.1 = 13 

(g) 

Rowland Bakke, 
Tmhof, Smith (1963) 
scintillators on 

March 1963 

9.2 

1.27 

.216 

340 km 

no maximum at 40 Mev 
J o ( 2 0) 0 n 

J 0 (40) ' 2 ’° 

J 0 (20) 

a 22 (Fig - 26) 


Atlas pods 






Jo (10) 

J 0 (60) 8 0 

Jo (10) 

j„ (60) “ 9 0 



Figure 69— The amplitude of the proton population 
change during the solar cycle for L = 1.4 and 1.6. 
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more) so that they can see proton-flux changes directly. From Table 2, experiments (c) 
and (d) are of this type. Although experiment (d) is not one experiment on one satellite, it 
is many identical experiments on essentially identical satellites and therefore falls in this 
category. 

2. Data on proton-flux changes obtained by comparing experiments on different satellites. 
Data of this type is given in experiments (a), (e) and (f) in Table 2. 

3. Data given by direct measurement of the proton energy spectrum in one experiment. Ex- 
periments (a), (b), (e), (f) and (g) are of this type. 

In comparing flux time change experiment of type 1 with the calculations, the agreement is 
good. Two experiments covering the period 1959 to 1962 give results consistent with the 
calculations. 

Comparing with type 2 data does not work as well. The observed changes of flux with time 
agree qualitatively with those predicted, but not quantitatively. The experimental flux ratios are 
smaller than the calculated ones. Two remarks are in order here: First, type 2 data usually has 
larger errors attached to it than type 1 because it involves the systematic errors of two experi- 
ments while type 1 data involves no systematic errors (only statistical errors are involved in the 
time variations in the flux). Because of this the experimental flux ratios using type 2 data are 
probably not as accurate as those using type 1 data. Secondly, all the comparisons of type 2 data 
involve Explorer IV data, so there is no independent check of experimental consistency. 

Of the five experiments of type 3 that help 
to understand the proton energy spectrum, the 
first two ((a) and (b) in I960) showed a modest 
sized peak at about E = 40 Mev. The more re- 
cent experiment (g) of Rowland et al (Reference 
24) shown in Figure 70 does not show such a 
peak. The solid curve on Figure 70 represents 
the normalized results of the present calcula- 
tions compared with the data of Rowland. The 
agreement is fairly good. Recent experiments 
((f) and (g)) show a large low-energy proton pop- 
ulation but cannot be compared with early ex- 
periments because the early ones did not go 
down to such low energies. Figures 41 through 
67 show that the calculated spectral shape should 
not change much during the period of these meas- 
urements, from 6 to 9 years solar-cycle time. 
The proton fluxes should increase but the shape 
of the spectrum stays nearly the same. Striking 
changes in spectral shape should occur on the 



ENERGY (Mev ) 

Figure 70— The energy spectrum of inner zone protons 
measured by Rowland et al, 1963 on an Atlas pad at 
L = 1.27, B = .216. 

Shown for comparison is the spectrum calculated ir> 
this paper for L = 1.27, B = .207, t ~ 9. 
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upswing of the solar cycle. The agreement with energy =25 Mev 

spectral shape in Figure 70 is not bad but the solar mm^* 

peaks in experiments (a) and (b) should not oc- x= 30 ° 

cur at t = 6. It may be that a neutron absorp- 
tion as suggested by Freden and White (Refer- 
ence 3) is required to explain the peak. 

All of the data presented so far are mirror 
point fluxes, M, (i.e.,the flux of protons mirror- 
ing per unit volume at one particular B and L). 

For low altitudes where the atmosphere changes 
rapidly with altitude, it is nearly correct to 
compare this flux with measured omnidirec- 
tional flux values, since, for these low altitudes, 
the omnidirectional flux is very nearly the same 
as the mirror-point flux. To show a more com- 
plete picture of the solar- cycle proton changes 
we have converted to omnidirectional fluxes 
using Equations 24 of Hess and Killeen (Refer- 
ence 27). Figure 71 shows a calculated R - \ 

map of the 25 Mev proton omnidirectional fluxes, J, at solar maximum and solar minimum. An 
isoflux contour is clearly at lower R for solar minimum. 



Figure 71 — An R - X. map of the omnidirectional flux 
of E = 25 Mev protons at solar max and solar min. 


CONCLUSIONS 

In conclusion, large changes in proton fluxes will take place during the solar cycle for those 
regions of b-L space corresponding roughly to minimum altitudes of about 300-800 km. Changes 
in proton spectra will occur also. Comparison of the calculations with available experimental in- 
formation are not conclusive. Some kinds of experimental data agree quantitatively with the cal- 
culations. Other data agrees qualitatively but not quantitatively. The crucial tests of the calcula- 
tions will involve measuring changes in the proton flux and spectrum as we approach solar maximum 
during 1966-68. We know experimentally (Reference 16) that the atmospheric density changes used 
in the calculation are reasonable. These calculated effects must take place with about the magni- 
tudes shown here unless there are features of the inner-belt protons which we do not now understand. 


(Manuscript received June 4, 1965) 
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